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This paper extends the results of a previous paper in which two models 
of a system of queues served in cyclic order were studied. One model is an 
exhaustive service model, in which the server waits on all customers in 
a queue before proceeding to the next queue in cyclic order. The other is a 
gating model, in which a gate closes behind the waiting units when the 
server arrives, and the server waits on only those customers in front of 
the gate, deferring service of later arrivals until the next cycle. 

In the present paper, the Laplace- Stieltjes transforms of the order-of- 
arrival waiting time distribution functions and, for the exhaustive service 
model, the mean waiting time for a unit arriving at a queue, are obtained. 

I. INTRODUCTION 

In a recent paper 1 we studied two models of a system of queues 
served in cyclic order: 

In each model, the ith queue is characterized by general service time 
distribution function //,-(•) and Poisson input with parameter \, . In 
the exhaustive service model, the server continues to serve a particular 
queue until for the first time there are no units in service or waiting in 
that queue; at this time the server advances to and immediately starts 
service on the next nonempty queue in the cyclic order. The gating 
model differs from the exhaustive service model in that when the server 
advances to a nonempty queue, a gate closes behind the waiting units. 
Only those units waiting in front of the gate are served during this 
cycle, with the service of subsequent arrivals deferred to the next cycle. 

In Ref. 1 we found, for the exhaustive service model, expressions for 
the mean number of units in a queue at the instant it starts service, the 
mean cycle time, and the Laplace-Stieltjes transform of the cycle time 
distribution function. 

In the present paper, we extend the analysis to obtain, for each 
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model, the Laplace-Stieltjes transform of the order-of-arrival waiting 
time distribution function and, for the exhaustive service model, the 
mean waiting time for a unit arriving at the ith queue. 

In Ref. 1 we defined a switch point as a time epoch at which the 
server finishes serving a queue; and we defined P, (fa , ■ • • , n N ) as the 
joint probability that at a switch point the server has just completed a 
visit at queue i (i = 0, 1, ■ ■ • , N) and Wi units are waiting in queue 
i + 1, n, units in queue i + 2, • • • , and n iV units in queue i + N. 

The central results of Ref. 1 were an iterative algorithm for the 
calculation of the probability generating functions 

0,(zi , • • • , Xn) = S • • • X) Pifai i • • * i n M l ••• xW 

n,-0 njv-0 

(t-0,1, ••• ,210 (1) 

and, for the exhaustive service model, an expression for the mean 
number «, of units waiting in queue i + 1 when the server completes a 
visit at queue i. In particular, it was shown for the exhaustive service 
model that these generating functions satisfy the functional equations 

0«(*i i • • • i *») = gi-i\Pi\1L x -+-( 1 - OJ i *x i • • ■ i *»-iJ 

+ y p<{l Xi-d - *-))^(o) - Pi-i<P, • • • , 0) 

(t = 0, 1, ••• ,iV) (2) 
and that n, is given by w, = m t /g { (l, • • ■ ,1), where 

^,. = ^±ip(0) p ~ Pf+1 (*-0,l f ••• ,N) (3) 

A 1 — p 

and where X,- is the rate of arrivals of units at queue i, p.- is the traffic 
intensity at queue i, 



it 



x = E v. 

• -0 

.V 

p = £ p.. 
P(0) = Ip,(o, ••• ,0), 



and /3,-(«) is the Laplace-Stieltjes transform of the distribution function 
of the length of the busy period at queue i. Equations (1), (2), and (3) 
appear in Ref. 1 as equations (3), (5), and (34), respectively. 
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The distribution generated by g i (x 1 , • • • , x N ) is denned with respect 
to a Markov chain imbedded at the switch points. An analysis by 
Takacs 2 of the exhaustive service model for the special case of two 
queues is based on a Markov chain imbedded at the set of service com- 
pletion points. Clearly, the set of switch points is a proper subset of 
the set of service completion points. Our use of switch points instead of 
service completion points enabled us in Ref . 1 to analyze the multiqueue 
model with about the same degree of mathematical complexity as 
Takacs required for the analysis of the 2-queue model. On the other 
hand, one would expect that our use of a chain imbedded in a "smaller" 
set of points would result in a corresponding loss of useful information. 

Takacs' analysis yielded waiting time results for the 2-queue model. 
At the time Ref. 1 was written, it was not apparent to us that our 
method of analysis provided enough information to enable us to obtain 
corresponding waiting time results for the multiqueue model. Accord- 
ingly, we concentrated on the cycle time, a quantity that seemingly 
gives the same kind of information as the waiting time. Unfortunately, 
we did not have complete freedom in choosing a precise definition of the 
quantity we would call cycle time. The mathematical formulation of 
the model dictated that the cycle time for queue i be denned, roughly 
speaking, as the length of time between two successive instants at which 
the server completes service at queue i, without regard to whether or 
not the server is continuously busy throughout this time interval. This 
definition introduces, among others, the following difficulty in the inter- 
pretation of realized values of the cycle time: A long cycle time could 
have resulted either from heavy congestion or from no congestion. 

No such ambiguities exist with respect to interpretation of the waiting 
time, which is simply the elapsed time from the arrival instant of a 
unit to the instant at which service on this unit begins. Therefore, we 
would like to obtain waiting time results. Furthermore, we would like 
to obtain these results, if possible, without directly extending the pre- 
vious analysis to include the entire set of service completion points. 

In the present paper, we obtain the desired waiting time results 
without recourse to a complicated reformulation of the original analysis 
based on the complete set of service completion points. Rather, to obtain 
the waiting times at queue i, we use the generating function 
g t -i(xi , • • • , x N ), calculated in Ref. 1, to append to the original set of 
switch points only those service completion points that correspond to 
departures from queue i; and this is sufficient for our waiting time 
calculations. 

The essence of the method is to calculate the probability generating 
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function of the number of units left in queue i by an arbitrary de- 
parture from queue i, using only the (known) probability generating 
function of the number of units waiting in queue i when the server 
arrives. The Laplace-Stieltjes transform of the order-of-arrival waiting 
time distribution function for units at queue i is then easily obtained 
by a standard argument. 

The preceding discussion refers mainly to the exhaustive service 
model, which was discussed in detail in Ref. 1. The gating model was 
shown to be characterized by equations that are essentially the same 
as those of the exhaustive service model, and was therefore not devel- 
oped in detail. In the present paper, the waiting times for the gating 
model will also be discussed. 

II. THE M/G/l QUEUE WITH SERVER VACATION TIMES 

In preparation for calculation of the waiting times in the exhaustive 
service model, we first consider the following generalization of the 
M/G/l queue: 

As usual, the server serves the queue continuously as long as there 
is at least one unit in the system (waiting or in service). When the 
server finishes serving a unit and finds the system empty, however, it 
goes away for a length of time called a vacation. At the end of the vaca- 
tion the server returns to the queue, and begins to serve those units, if 
any, that have arrived during the vacation. If the server finds the 
system empty at the end of a vacation, it immediately takes another 
vacation, and continues in this manner until it finds at least one waiting 
unit upon return from a vacation. 

Let X k (k = 1, 2, • • •) be the number of units left behind by the /bth 
departing unit. Then 

P[X k+1 = n) = "JtPlX, = v}P{X k ^ =n\X k = v\ 

(fc = 1,2, ... ;n = 0,1, •••)• (4) 

Let P(j) be the probability that at the end of a vacation the server 
finds j ^ units waiting for service. If the arrival rate and the service 
time distribution function are denoted by X and H(-), respectively, then 

P{X k+l - n | X k - v > 0} - / o " (w ( ^ ) " 1 + - t ;)! exp ( ~ X * } dm) 

(n ^ v - 1) (5) 
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and 
P[X M =n |X. = 0} = %T*% /; 5r ^J I exp(-X0 

<*#(*). (6) 

The expression P(j)/[1 — P(0)] in equation (6) is the conditional prob- 
ability that when the server starts serving the queue there are ; units 
waiting, given that at least one unit is waiting. 

When the traffic intensity p is less than unity (p = \h, where h is 
the mean service time), there exists a unique distribution 

7r„ = UmP{X k =n} (n = 0, 1, •••) (7) 

*-«00 

that satisfies both the normalization equation 

£».-! (8) 

n-0 

and the limiting set of equations obtained from equation (4), 

+ 2 r. {" (n+^r-l)! 6XP ( ~ X ^ dffft) (n = °' l| ' ' ' } - (9) 
Define the probability generating functions 

M - Et/ (io) 

n-0 

and 

*go - £pov. (id 

Substitution of equation (9) into equation (10) yields, after some 
manipulation, 

where r?(-) is the Laplace-Stieltjes transform of the service time distri- 
bution function H(-). Observe that the expression [\fr(z) — P(0)]/[1 — 
P(0)] is the probability generating function of the number of units 
waiting when service commences. 
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The unknown probability t is determined from equation (12) by the 
normalization condition /(l) = 1. Application of l'Hospital's rule to 
equation (12) yields 

*• "Too ( " p) - (13) 

Thus, the probability generating function /( • ) of the number of units left 
behind by an arbitrary departure, and the probability generating func- 
tion \p(-) of the number of units waiting at the end of a vacation are 
related as follows: 

KX) ~ x- v (X- Xx) W) ' (14) 

We now apply a standard argument to obtain the Laplace-Stieltjes 
transform co(-) of the order-of-arrival waiting time distribution function 
from the generating function (14). 

Let F(-) be the distribution function of an arbitrary unit's sojourn 
time, defined as the elapsed time between the unit's arrival and de- 
parture epochs, and denote by <£(•) the Laplace-Stieltjes transform 
of F(-). Since the sojourn time is the sum of the waiting time and the 
service time, and since these latter times are independent, therefore 

*(«) = co(s)r?(s). (15) 

When units are served in their arrival order, each departing unit 
leaves behind it precisely those units that arrived during the departure's 
sojourn time. Further, these remaining units arrived according to a 
Poisson process, independent of the sojourn time. Therefore, the prob- 
ability t„ that a departing unit leaves behind n other units is 

*• - f m ^r ex P (-*0 dF ®- < 16 > 

Jo »*' 

Substitution of equation (16) into (10) gives the well known and funda- 
mental relation 

m = ^(x - \x). (17) 

Equations (14), (15), and (17) together give the Laplace-Stieltjes 
transform co(-) of the waiting time distribution function in terms of the 
probability generating function \f/(-) of the number of units waiting 
at the end of a vacation: 
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Note that for the ordinary M/G/l queue, in which the vacation ends 
immediately whenever a unit arrives and finds the server idle, ip(x) = x 
and equation (18) reduces to the well known Pollaczek-Khinchin 
formula, 

«•« — *llw < (19) 

s — X + Xi?(s) 

Finally, the mean wait for service W = — w'(0), obtained from 
equation (18), is given by 

W = 20^-p) + 2X7(1)' (20) 

The first term on the right side of equation (20) is identical with the 
mean waiting time in the ordinary M/G/l queue, as would be obtained 
directly from the Pollaczek-Khinchin formula (19). The second term 
in equation (20) represents the component of the mean wait that arises 
because of the variability in the number of units waiting when service 
begins. Although service in arrival order was assumed in its derivation, 
equation (20) is valid for any order of service that is independent of the 
service times. 

Note that the result (14) is true regardless of any relationship be- 
tween the vacation lengths and the arrival process, whereas equations 
(18) and (20) are valid only when the vacation lengths are determined 
without regard to the arrival process. For example, if one were to 
consider a mechanism such that service begins as soon as a fixed number 
j (j ^ 2) units are waiting, then equation (14) with \[/(x) = x j would 
correctly give the probability generating function of the number of units 
in the system just after a service completion epoch. On the other hand, 
equations (18) and (20) would not apply, because (16), and therefore 
(17), would no longer be true. (For if ${x) = x', then the first departing 
unit would always leave behind at least those j — 1 other units that 
were present when service commenced. Thus, for the first departing 
unit, t b = for n < j — 1, and this contradicts the assumption (16) 
if j ^ 2.) 

III. LAPLACE-STIELTJES TRANSFORM OF WAITING TIME DISTRIBUTION 
FUNCTION FOR EXHAUSTIVE SERVICE MODEL 

We now proceed to apply the results of Section II to the analysis of 
waiting times in the exhaustive service model. In essence, the "vaca- 
tion time" of Section II is the length of time that the server spends idle 
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or working on other queues before registering a switch point at queue 
i — 1 and beginning service on queue i. 

For want of a better word, let us define as a supercycle the elapsed 
time between the arrival epoch of a unit at any queue when the system 
is completely empty, and the first instant at which the whole system 
again becomes empty. Then units that arrive at queue i can be classified 
into two exclusive and exhaustive categories: 

(1) arrivals at queue i that either initiate a supercycle or occur 
during the 1-busy period generated by an arrival at queue i that 
initiated a supercycle; or 

(2) all other arrivals at queue i. 

Equivalently, units in category (1) are those arrivals at queue i whose 
service begins prior to the occurrence of the first switch point of a 
supercycle, whereas units in category (2) are arrivals at queue * whose 
service times begin after the first switch point of a supercycle. 

Consider now the waiting times of units that arrive at queue i. Those 
units in category (1) are served during a busy period originated by one 
unit. Therefore, the Laplace-Stieltjes transform w" } (-) of the order- 
of-arrival waiting time distribution function for units at queue i that 
belong to category (1) is given by the Pollaczek-Khinchin formula (19) : 

«!"(.) = S( ' 1 ?> , , (21) 

s — X, • + \iiji(s) 

where p, , X,- , and Vi(') are the corresponding quantities in equation (19) 
defined now with respect to queue *. 

Units in category (2) are served during a busy period originated by 
those units waiting in queue i when the server leaves queue % — 1. 
Let ^i(-) be the probability generating function of the number of units 
waiting in queue i when the server leaves queue i — 1; ^.(') is the 
probability generating function of the number of units waiting for 
service in queue i when the server finishes a vacation, and is given by 

*« - g:!'!: ::::!) ■ < 22 > 

[Note that ^'(1) = »<_i = fhi-i/gt-x (1, 1, • • • , 1).] Thus, the Laplace- 
Stieltjes transform coJ 2) (-) of the waiting time distribution function for 
units in category (2) is given by equation (18) : 



"fw-ifet 1 -*^ 



1 " (23) 



_ s — X< + \iVi(s) 
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Let p< 1} be the proportion of all arrivals at queue i that are in cate- 
gory (1). The mean number of units that arrive at queue i during an 
interval of length t is \ t t. The probability that an arbitrary arrival at 
queue i finds the whole system empty is 1 — p (p = p + • • • + Pn), 
so that \it(l — p) is the mean number of arrivals at queue i that initiate 
a supercycle during any elapsed time t. The mean number of units 
served at queue i during the 1-busy period generated by each such 
arrival is (1 — p,) _1 , and hence the mean number of units in category 
(1) served at queue i during an elapsed time t is \ { t(l — p)/(l — Pi)- 
Therefore, the probability is [X,-<(1 - p)/(l - p,-)]/M that an arbitrary 
arrival at queue i is in category (1); that is, 

o> = 1^_P (24) 

1 — Pi 

and the probability p- 2) = 1 — pj" that an arbitrary arrival at queue i 
is in category (2) is 

« = P_ZL-2±. (25) 

1 — Pi 

The Laplace-Stieltjes transform co.(-) of the waiting time distribution 
function for an arbitrary unit at queue i is the weighted sum of the 
transforms for each category: 

Ui (s) =pl. 1) w » ) ( S )+ ? >! 2> a ) , (2, (s). (26) 

Finally, equation (26) becomes, with the help of equations (21) through 
(25) and equation (3), 

M 1 ~ P 

s — X, • + \iVi(s) 



U-i(l, 1, 



i-A-L .« i ... n-.,.^ 



i) - 0.-. 



+ s 



\P(o) L y,_iv ' ' ' y,_1 v h ' 

(i = 0,l, ••• ,iV). (27) 

Inherent in equation (27) is the assumption that units in queue i are 
served in their arrival order, but no assumption is made regarding the 
order of service of units in other queues. If at each queue units are 
served in their arrival order, then the waiting time distribution function 
for an arbitrary unit, without regard to the identity of the queue in 
which it is served, has Laplace-Stieltjes transform w(-) given by 

w(s) = X" 1 i>,«j,(s). (28) 
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IV. MEAN WAITING TIMES FOR EXHAUSTIVE SERVICE MODEL 

Denote by W t the mean wait for service suffered by units arriving 
at queue i. The mean wait for service for units in category (1) is 
[X,-77$'(0)/2(l — pi)]; the mean wait for service for units in category (2) is, 
in analogy with equation (20), [X,i7{'(0)/2(1 - p,)] + [#'(1)/2X<#(1)]. 
Weighting these values according to equations (24) and (25), respec- 
tively, we have 



Wi 2(1 - p.-) + tf(l) 2X,(1 - p,) 
In equation (26) of Ref. 1 we defined 



(* = 0, 1, 



N). (29) 



w,(fc) = — gfa , 



> x n) 



a - o, i, 

and m, = m,(l). Let us also define 
X(l - p) d" 



N; k = 1, 



N) (30) 



rriiij, k) = 



gifri , ••- ,x N ) 



P(0) dx,dx k *, « w -i 

(i = 0,l, ■-• ,N; j- 1, •■■ ,N; k = 1, 
Then it follows from equation (22) and these definitions that 

#'(1) 



*'(!) 



P(0) m.-^l, 1) 
X(l - p) ^,-,(1) 



V)- (31) 
(32) 



Using equations (3) and (32), we can rewrite (29): 

(* = 0, 1, 



W{ W'(0) , ™,-:(l,D 



2(1 - p.) 2X?(1 - p.) 



,A0- (33) 



It remains to calculate the quantity ra,_i (1, 1) in (33). To this end, we 
define 



J««-£/»*(EWi-*o)| 

(t-0,1, •■ 

and 

(t = 0,l, •■- ,N; j =1, 



•-XN-1 

N; k = \, 



-XAT-l 

,N; fc-1, 



, #) (34) 



,A0- (35) 
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Note that in terms of the given parameters, 



1 — pi 



and 



M, k) = \. + ,A, +t {l v '[_ (0) p)a (37) 

where h { is the mean and ??,•(•) the Laplace-Stieltjes transform of the 
service time distribution function for a unit at queue i. 

We proceed to calculate w,(l, 1) in the same way we calculated 
ra,(l) in Section VII of Ref. 1. Differentiating twice through equation (2) 
and setting Xi = ■ ■ ■ = x N = l,\ve obtain the three-dimensional set of 
linear equations 

m«0', *0 - x( p ( ~ } p) M<-iQ)M, fc) + (i - p)xAG, fc) 

+ A(j)ACb)A l . x ai i) + (i - *(# - j))j8^)^-i(i, ; + i) 

+ (i - s(at - fc))A0)m,- l a, * + i) 

+ (i - «(JV - i))(i - «(iv - fc))«,-,0 + 1, k + 1) 

(< - 0, 1, ••• ,JV; j-1, -•• ,N; k = l,---,N) (38) 

where 5 (x) = 1 if a; = and 8 (x) = if x ^ 0. Using equation (3) and 
combining the first two terms on the right side of equation (38), we can 
write 

ffiiij, k) = \,(1 - p.)/3,(i, k) + /3,(;)/3.(/c)w,- 1 (l, 1) 

+ &(k)mi-i(X, i + 1) + /3.0')w,-i(l, fc + 1) + w<-i(j + 1, fc + 1) 
(t - 0, 1, • • ■ , N; j = l,---,N; k = 1, • • • , N) (39) 

where all undefined terms are taken to be zero. (The functions ffl t (j, k) 
are defined only for J, k = 1, • • • , N.) It is required to solve this set of 
%N(N + l) 2 independent linear equations for the m,(l, 1). [Note that 
ffiiij, k) = fhi(k, j).] 

Successive substitution into the last term on the right side of equa- 
tion (39) gives 

w,(l, k) = E *<-'(! - *->)&-»(! + ", fc + ") 

+ E A-.(l + ")/3,-,(fc + r)«i-»-p(li 1) 
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+ £&-r(fc+")^i-l-r(l,2+ J ') 

v = 

+ £&_,(1 + ,)m,_ 1 _,(l,/c + 1 +v) 

? = 

(t = 0, 1, ••• ,AT; /c = 1, •■• ,N) (40) 

where each sum is continued as long as the terms are denned. 

The set (40) consists of N(N + 1) independent linear equations. 
Unfortunately, it does not appear that further algebraic simplification 
is likely. However, for particular values of the parameters and reason- 
able values of N, numerical solution should not be difficult. 

Therefore, to calculate the mean wait W ,• for service at queue 
i (i = 0, 1, • ■ ■ , N) for any particular values of the basic parameters, 
simply solve the set (40) of N(N + 1) linear equations numerically, 
and use the resulting value of rn t -i (1, 1) in equation (33). Note that 
these calculations for the mean waiting times require no iteration, since 
neither generating functions nor state probabilities appear. This last 
observation is remarkable in light of the complicated iteration process 
(discussed in Ref. 1) underlying the derivation of these results. Thus, 
despite the complicated derivation, calculations do not seem 
impractical. 

In the particular case of two queues (iV = 1), only N(N + 1) = 2 
simultaneous equations must be solved to find ra<(l, 1), and an algebraic 
solution is easily obtained. For N = 1 equation (40) gives 

_ „ n , 2 A.-^^OXl ~ P,) 2 + X.^'(0)p*-i 
m,_,(l, 1) = X,- (1 _ p _ i)2(1 _ p)2 _ p2 _ )p2 

and hence, for two queues, 

^ w(Q) x,-. 1 i ? ;: 1 (0)(i - p,) 2 + hv7(Q) P U 

Wi ~ 2(1 - Pl ) + 2(1 - p,)[(l - P.-i) 2 (l - P.) 2 " P 2 -iP 2 ] 

(i = 0, 1). (41) 

Our result (41) is in agreement with previous results of Takacs, 2 Avi- 
Itzhak, Maxwell and Miller, 3 and Eisenberg. 4 

Although service in order of arrival was assumed throughout, the 
results for the mean waiting time are valid for any order of service that 
is independent of the service times. 

V. LAPLACE-STIELTJES TRANSFORM OF WAITING TIME DISTRIBUTION 
FUNCTION FOR GATING MODEL 

Turning to the gating model, we now study briefly the distribution 
of waiting times for units served in order of arrival at the ith. queue. 






CYCLIC QUEUES 411 

As with the exhaustive service model, we first calculate the probability 
generating function of the number of units left in queue * by an arbi- 
trary unit departing from queue i; and using the same arguments, we 
obtain from this generating function the Laplace-Stieltjes transform 
of the waiting time distribution function. As in Ref. 1, the notation 
for the gating model is the same as for the exhaustive service model; 
<7,_i (xi , • • • , Xat) and related probabilities are defined and calculated as 
described in Section IX of Ref. 1. 

Let Ti(j) be the conditional probability that an arbitrary departure 
from queue i leaves behind it j units in queue i, given that this de- 
parture did not arrive when the system was completely empty. Then 

*» - tt^Sml SW*™ *<-*-« 

(i = 0, 1, •••) (42) 

where P,_i(w)/[1 — P { - x {0)} is the conditional probability that n ^ 1 
units are waiting in queue i when the gate closes, given that at least 
one unit is waiting; and \/n is the probability that a departing unit is 
fcth in line for service (k = 1,2, • • • , n) given that n units are present 
at the closing of the gate. The integrand in equation (42) is taken to be 
zero when n — k > j. 

Following the argument for the exhaustive service model we see that 
the (conditional) probability generating function of the number of 
arrivals at queue i that occur during the waiting time of a departing 
unit (given that the departing unit did not find the system empty on 
arrival) is 

E7r,(i).T y A.(X.--X,T). 

1=0 

A simple calculation gives 

CO 

-9? — ttt - £ i p ( L - E arvrcx, - m). (43) 

i?<(X< — A,.r) tTi 1 — ".--i(0) n k -i 

The probability that an arbitrary departing unit did not find the system 
completely empty on arrival is p = S'T=o Kh t . Thus, after summing 
the geometric series in equation (43), the unconditional order-of-arrival 
waiting time distribution function for units served at queue i has La- 
place-Stieltjes transform w.(-) given by 
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M n v , y PUn) i [\<v<(8)] n - (x« - s) n 

«,(«) = (1 - P ) + P Zjj.p^Qnxr 1 a - X.- + M,(*) 

(i = 0, 1, ••• ,N - 1). (44) 

Equation (44) allows numerical calculation (and hence numerical 
inversion) of the transform «,(•)• Unfortunately, this procedure requires 
knowledge of the distribution {P < _ I (n)} l which is specified only through 
its generating function <7,_i (x, 1, • • • ,1). Thus, to obtain numerical 
results for the gating model one must solve two distinct problems in nu- 
merical analysis, numerical calculation of the {P t -_ x (n)} and then nu- 
merical inversion of the transform. Note that the first of these numerical 
calculations is not required for the exhaustive service model. The subject 
of numerical inversion of Laplace-Stieltjes transforms and probability 
generating functions (the latter being, in fact, a special case of the 
former) is important for the reduction to practice of these cyclic queuing 
models. However, it is a subject best treated separately, without regard 
to the particular applications at hand, and will not be discussed further 
here. 

VI. SUMMARY AND PROPOSALS FOR FUTURE WORK 

We have extended our previous study of cyclic queues to obtain 
waiting time results. In particular we have obtained, for both the ex- 
haustive service model and the gating model, the Laplace-Stieltjes 
transform of the waiting time distribution function for units arriving 
at the tth queue, when units at that queue are served in order of arrival. 
These transforms are given by equation (27) for the exhaustive service 
model and equation (44) for the gating model. Also, we have obtained 
for the exhaustive service model a formula (33) for the mean waiting 
time for units arriving at the ith. queue. Use of equation (33) requires 
calculation of the value M t -i (1, 1), which can be obtained in any 
particular case by numerical solution of the N(N + 1) linear equations 
(40). It is noteworthy that the calculation of the mean waiting time 
requires no iteration. 

The techniques used in this and our previous study might be useful 
in the analyses of priority queuing models and other cyclic queuing 
models that have important practical applications. Examples of the 
latter are: extensions of the present models to include arbitrary switching 
times and/or set up times; systems of queues served in arbitrary periodic 
order (of which cyclic order is a special case); and within-queue disci- 
plines other than service in order of arrival, such as service in random 
order. 
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